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We calculate the temperature dependence of the fluorescence Stokes shift and the fluorescence 
decay time in linear Frenkel exciton systems resulting from the thermal redistribution of exciton 
population over the band states. The following factors, relevant to common experimental conditions, 
are accounted for in our kinetic model: (weak) localization of the exciton states by static disorder, 
coupling of the localized excitons to vibrations in the host medium, a possible non-equilibrium of 
the subsystem of localized Frenkel excitons on the time scale of the emission process, and different 
excitation conditions (resonant or non resonant). A Pauli master equation, with microscopically 
calculated transition rates, is used to describe the redistribution of the exciton population over 
the manifold of localized exciton states. We find a counterintuitive non-monotonic temperature 
dependence of the Stokes shift. In addition, we show that depending on experimental conditions, 
the observed fluorescence decay time may be determined by vibration-induced intra-band relaxation, 
rather than radiative relaxation to the ground state. The model considered has relevance to a wide 
variety of materials, such as linear molecular aggregates, conjugated polymers, and polysilanes. 

PACS numbers: 



I. INTRODUCTION 

The dynamics of excitons in chain-like systems, like 
conjugated polymers^ polysilanes, 2 and molecular J- 
aggregatesjMi have attracted much attention over the 
past decades. The dynamics in these systems result from 
the complicated interplay between various processes, in 
particular from relaxation to the ground state (popula- 
tion relaxation), from migration of the exciton to sites 
with different excitation energy (population redistribu- 
tion), and from relaxation caused by nuclear displace- 
ments (exciton self-trapping). Population relaxation is 
distinguished in radiative and nonradiatiave channels. 
With regards to migration, two limiting situations are 
usually considered. The first one is the case of incoher- 
ent energy transfer (Forster transfer), which takes place 
between strongly localized excitations on the chain. This 
is also referred to as hopping transfer in a disordered site 
energy distribution. In the other limit, one deals with 
weakly localized exciton states, which may extend over 
many repeat units of the chain. Transitions between such 
band-like exciton states are possible due to their scat- 
tering on lattice vibrations and are often referred to as 
intraband relaxation. 

The net effect of the above processes may be probed 
by various optical techniques, of which steady-state and 
time-resolved fluorescence spectroscopy are the most fre- 
quently used ones. Common quantities extracted from 
such experiments are the decay time of the total flu- 
orescence intensity following pulsed excitation and the 
steady-state as well as the dynamic Stokes shift. As most 
of the dynamic processes mentioned above are influenced 



by temperature in a different way, one often also probes 
the temperature dependence of the fluorescence. 

The temperature dependence of the fluorescence 
decay time has drawn particular attention in the case 
of molecular J-aggregates of polymethine dyes. In 
these systems, the Frenkel exciton states are delocalized 
over tens of molecules (weakly localized excitons). 
The coherent nature of the excitation extended over 
many molecules leads to states with giant oscillator 
strengths, which scale like the number of molecules over 
which the exciton state is delocalized. These so-called 
superradiant states lie near the bottom of the bare 
exciton band, and, especially at low temperatures, 
lead to ultrafast spontaneous emission (10's to 100's of 
ps) .SiSigiiLiSiiiiiSiiiili Upon increasing the tempera- 
ture, these systems typically exhibit an increase of the 
fluorescence lifetime, which is a trend that is highly 
unusual for single-molecule excitations and is intimately 
related to the extended nature of the exciton states. 
This unusual temperature dependence was first observed 
for the J-aggregates of pseudoisocyanine (PIC)f£i2*i2i 
and later has been confirmed for other types of 
J-aggregates, in particular, 5,5',6,6'-tetrachloro-l,l'- 
diethyl-3,3'-di (4-sulfobutyl)-benzimidazolocarbocyanine 
(TDBC)ri l,l'-diethyl-3,3'-bis(sulfopropyl)-5,5',6,6'- 
tetrachlorobenzimidacarbocyanine (BIC)^ and 3,3'- 
bis(sulfopropyl)-5,5'-dichloro-9-ethylthiacarbocyanine 
(THIATS)iii^ The slowing down of the aggregate 
radiative dynamics is usually attributed to the thermal 
population of the higher-energy exciton states, which in 
J-aggregates have oscillator strengths small compared to 
those of the superradiant exciton states£& In spite of the 
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basic understanding that the redistribution of the initial 
exciton population over the band states plays a crucial 
role in this problem, the theoretical efforts to describe 
this redistribution and to fit all details of this behavior 
have not been fully successful so f aB jL±iii§iiLi& (see the 
discussions presented in Ref. ITflh . In particular, as we 
have found from our previous study on homogeneous 
model aggregates^ the initial excitation conditions 
seem to play a crucial role in the interpretation of the 
measured fluorescence lifetime: under certain conditions 
this lifetime probes the intraband relaxation time scale, 
rather than the superradiant emission time. 

The Stokes shift of the fluorescence spectrum in lin- 
ear exciton systems has been considered previously by 
various authors; studies of its temperature dependence 
are rare. One of the first studies concerned the Stokes 
shift in polysilanes at liquid-helium temperatures, which 
was modelled in a phcnomcnological way by assuming 
that intraband relaxation is determined by one energy- 
independent relaxation rate in combination with the 
number of available lower-energy exciton states in a dis- 
ordered chain^ Relaxation by migration to the lowest- 
exciton state available on a linear chain also forms the 
main ingredient of the theoretical study of Chernyak et 
al^H on the relation between the fluorescence lineshape 
and the superradiant emission rate deep in the red wing 
of the density of states of disordered J-aggregates at cryo- 
genic temperatures. For 7r-conjugated polymers the time- 
dependent shift of the luminescence spectrum (dynamic 
Stokes shift) has been modeled by assuming excitons lo- 
calized on a few repeat units, which migrate spatially as 
well as energetically through Forster transfer £±22^ It 
was concluded that the Stokes shift in these systems can- 
not be explained from nuclear displacements, as would 
be the case for single molecules, and that the migration 
process plays a crucial rolei^I Also for J-aggregates, with 
strongly delocalized exciton states, the believe is that the 
Stokes shift induced by nuclear displacements is small, in 
fact much smaller than the Stokes shift of their single- 
molecule constituents. The explanation lies in the fact 
that for a delocalized excitation the weight of the excita- 
tion on each molecule of the chain is small, which leads to 
a small nuclear displacement on each molecule. The mea- 
sured Stokes shift differs for various types of J-aggregates, 
for instance for PIC a shift can hardly be detected^ while 
for TDBCjii BIC,±2i and THIATSi^ a clear shift can be 
observed. As far as we are aware, the temperature depen- 
dence of the Stokes shift was only measured for THIATS 
aggregates^ It shows an interesting non-monotonic be- 
havior, analogous to the one found in disordered narrow 
quantum wells,— where it finds its origin in thermally 
activated escape from local minima in the random po- 
tential. 

The goal of the present paper is to model the tem- 
perature dependence of the dynamics of weakly local- 
ized excitons in linear chains and to establish the ef- 
fect on the fluorescence lifetime and Stokes shift. We 
will be mostly interested in temperatures up to about 



100 K, where scattering on acoustic phonons is the 
dominant scattering mechanism. The current paper is 
an extension of our previous work, where we studied 
the temperature dependence of the fluorescence lifetime 
and restricted ourselves to ideal (homogeneous) Frenkel 
chains^ Here, we consider a more detailed model, which 
includes on-site (diagonal) disorder. It is well-known that 
this model provides a good basis for understanding the 
complex (linear and nonlinear) optical dynamics in J- 
aggregatesi 7 i 8 i 25 i 26 i 27 i 28 i 29 We will take into account the 
following factors that seem to be essential under common 
experimental conditions: (i) localization of the exciton 
states by the site disorder, (ii) coupling of the localized 
excitons to the host vibrations (not only to the vibrations 
of the aggregate itself as was done in Refs. ^| and fiH) . 
(iii) a possible non- equilibrium of the subsystem of local- 
ized Frenkel excitons on the time scale of the emission 
process, and (iv) the nature of the excitation conditions 
(resonant versus non-resonant). We use a Pauli master 
equation to describe the evolution of the populations of 
the localized exciton states and the intraband redistri- 
bution of population after the initial excitation. Previ- 
ously, such a master equation was also used to model 
the exciton dynamics in disordered quantum wells 2 ^ and 
polysilane films^ As observables, we focus on the Stokes 
shift of the fluorescence spectrum and the decay times of 
the total as well as the energy dependent fluorescence 
intensity. 

The outline of this paper is as follows. In Sec. |H] we 
present the model Hamiltonian of Frenkel excitons with 
diagonal disorder, the main effect of which is localiza- 
tion of the exciton states on finite segments of the aggre- 
gate. We briefly reiterate the basic facts concerning the 
structure of the exciton cigenenergies and eigenfunctions 
close to the band bottom, which is the spectral range 
that mainly determines the exciton optical response and 
dynamics. The Pauli master equation that describes the 
transfer of populations between the various exciton eigen- 
states, is introduced in Sec. IIIII The numerical solution 
of this equation under various conditions is obtained and 
used in Sec. El to study the steady-state fluorescence 
spectra and the temperature dependence of the Stokes 
shift and in Sec. to study the temperature dependent 
fluorescence lifetime. The results are also discussed in 
terms of back-of-the envelope estimates based on the low- 
energy exciton structure. Finally, we summarize and con- 
clude in Sec. IVII 



II. DISORDERED FRENKEL EXCITON MODEL 

We consider a generic one-dimensional Frenkel exciton 
model, consisting of a regular chain of N optically active 
sites, which are modeled as two-level systems with par- 
allel transition dipoles. The corresponding Hamiltonian 
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reads 

N N 

H = Y^ tn\n){n\ + 4m \n)(m\ , (1) 

n— 1 n,m 

where \n) denotes the state in which the nth site is ex- 
cited and all the other sites are in the ground state. The 
excitation energy of site n is denoted e„. We will account 
for energetic disorder by assuming that each e„ is taken 
randomly, and uncorrelated from the other site energies, 
from a gaussian distribution with mean e and standard 
deviation a. Hereafter, eo is set to zero. The hopping 
integrals J nm are considered to be nonrandom, and are 
assumed to be of dipolar origin: J nrn = — J/\n — m\ 3 
(J nn = 0). Here, the parameter J represents the nearest- 
neighbor coupling, which is positive for the systems of 
our interest, namely those that have the optically domi- 
nant states at the bottom of the exciton band. Molecular 
J-aggregates are prototype examples of such systems. Di- 
agonalizing the N x N matrix H nm — (n\H\m) yields the 
exciton eigenenergies and wavefunctions. In particular, 
the vth eigenvalue E v (with v = 1 , . . . , N) is the energy 
of the exciton state \v) = Yln=i Vtml 71 ); where ip un is the 
nth component of the J^th eigenvector. 

As has been shown in Ref.l3ll in the absence of disorder 
(er = 0, the case of a homogeneous chain) the eigenvectors 
in the presence of long-range dipolar interactions with an 
accuracy of the order of TV -1 agree with those for the case 
of nearest-neighbor hopping: 

/ 2 \ 1/2 N 
|fc)= (ivTl) £sin(ifn)|n) , ( 2 ) 

^ ' n— 1 

where we introduced the wavenumber K — 7rk/(N + 1) 
and k = 1, . . . , N is used as quantum label for this ho- 
mogeneous case. These states are extended (delocalized) 
over the entire chain. It turns out that the state fc = 1 is 
the lowest (bottom) state of the exciton band. Close the 
bottom (k <C N, i.e., K <C 1) and in the limit of large 
N, the exciton dispersion relation reads^i 

E k = -2A04J+J[^-lnK S j K 2 . (3) 

Furthermore, assuming that the chain is short compared 
to an optical wavelength, the oscillator strengths of the 
states \k) close to the bottom of the band are given by 

F fc = -g-ff sinirnV^-^ * . (4) 
N +1 \ ^ J N + 1 K 2 y J 

Here, the oscillator strength of a single molecule is set 
to unity. According to Eq. J3J, the lowest state k = 1 
(with the energy E\ = —2.404 J) accumulates almost the 
entire oscillator strength, Fx — 0.81 (A^+ 1). Its radiative 
rate is thus given by is 71 = 70^1 = 0.8l7o(A^ + 1), 
i.e., roughly N times larger than the radiative rate 70 



of a monomeriS The k = 1 state is therefore referred 
to as the superradiant state. The oscillator strengths 
of the other odd states (k — 3,5, ...) are much smaller, 
Fk = Fi/k 2 , while the even states (k = 2,4,...) carry 
no oscillator strength at all, Fk = 0. As a consequence, 
the exciton absorption band occurs at the bottom of the 
exciton band, red-shifted with respect to the monomer 
absorption band. This is characteristic for J-aggregates. 




FIG. 1: (a) Density of states p(E), (b) oscillator strength 
per state F(E), and (c) absorption spectrum A(E) calculated 
for different disorder strengths: a — 0.1J (solid), a = 0.3 J 
(dashed), and a = 0.5 J (dotted). The disorder results in a 
tail of states below the bare exciton band edge Ei — — 2.404J, 
which is more pronounced for larger degree of disorder. Each 
spectrum for the oscillator strength per state has a well de- 
fined maximum and shows that the states in the tail carry 
most of the oscillator strength. This spectrum also tends to 
widen upon increasing a /J. The absorption spectra simply 
reflect the fact that A{E) = p{E)F{E). 

In the presence of disorder [a ^ 0), the exciton wave 
functions become localized on segments that are smaller 
than the chain length N. One of the important conse- 
quences of this localization is the appearance of states be- 
low the bare exciton band bottom E\ = — 2.404J; these 
states form a tail of the density of states and in fact 
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carry most of the oscillator strength. As a consequence, 
the linear absorption spectrum of the exciton system is 
spectrally located at this tail. All these properties are il- 
lustrated in Fig. [I] where the density of states, p(E), the 
absorption spectrum, A(E), and the oscillator strength 
per state, F{E) — A(E)/p(E), are depicted for three val- 
ues of the disorder strength: a = 0.1J, 0.3J, and 0.5J. 
These quantities have been calculated in the standard 
way using numerical simulations and the definitions: 

p( E ) = ^ (y, s ( e - E »)^ > ( 5a ) 



= iv \£ (J>-J S ( E - E »)) > ( 5b ) 

where the angular brackets denote the average over the 
disorder realizations. The statistics was improved using 
the smoothening technique developed in Ref. |U 

Despite the fact that the tail of the density of states 
does not show any spectral structure, it has been shown 
that the exciton wave functions and energy levels in 
this spectral region do exhibit a specific (local) struc- 
turei 17 i 31 ' 34 i 35 i 36 This structure is revealed by plotting 
the wave functions obtained for a particular realization 
of the disorder (see Fig.[3]for the case of a = 0.1J). The 
tail states (filled in black) have an appreciable amplitude 
only on localized segments of a typical size N* (local- 
ization length; in the current example N* w 50). Some 
of them have no nodes within the localization segments 
(states 1-6 and 8 in Fig.0 with the states counted start- 
ing from the lowest one). Such states can be interpreted 
as local excitonic ground states. They carry large oscil- 
lator strengths, approximately N* times larger than that 
of a monomer, and thus mainly contribute to the exci- 
tonic absorption and emission. The typical spontaneous 
emission rate of these states is 7^ joN*. 

Some of the local ground states have partners local- 
ized on the same segment; examples are the doublets of 
states (2,9) and (6,10) in Fig. [21 These partner states 
have a well defined node within the localization segment 
and can be assigned to the first (local) excited state of 
the segment. Their oscillator strengths are typically sev- 
eral times smaller than those of the local ground states. 
Sometimes, but less frequently, a local manifold contains 
three states, such as the triplet (1,7,11), with the third 
state being similar to the second excited state of the seg- 
ment and having an oscillator strength small compared 
to that of the local ground state as well. The rest of 
the local ground states (see the states 3, 4, 5 and 8) do 
not have well defined partners, because the latter (higher 
in energy) are extended over a few (adjacent) AT*-site 
segments (the states 12, 13, 14 and 15). The oscillator 
strengths of these high-energy state are also small com- 
pared to those of the local ground states. Such states 
form mixed manifolds. 




Site number 

FIG. 2: Exciton wave functions tp vn and energy levels E v in 
the vicinity of the bottom of the exciton band for a particular 
realization of the disorder at a = 0.1J. The states are ob- 
tained by numerically diagonalizing the exciton Hamiltonian 
H nm for a chain of 500 sites. The baseline of each state repre- 
sents its energy in units of J. The wave function amplitudes 
are in arbitrary units. It is seen that the lower states (filled 
in black) are localized on segments of the chain with a typ- 
ical size small compared to the chain length. Some of these 
localized states can be grouped into local manifolds of two 
or sometimes three states that overlap well with each other 
and overlap much weaker with the states of other manifolds 
[see the doublets of states (2,9) and (6,10), and the triplet 
(1,7,11) (states counted starting from the lowest one)]. The 
higher states (filled in gray) are more extended and cover sev- 
eral segments at which the lower states are localized. 



The mean energy spacing between the levels of a seg- 
ment represents the natural energy scale in the tail of 
the density of states. The local spectral structure corre- 
sponding to this spacing, however, turns out to be hidden 
in the total density of states, because the mean absolute 
energy difference between the local ground states of dif- 
ferent segments is approximately 1.5 times as large as the 
energy spacing between states within a single segment M> 
In spite of the fact that this local structure is not visible 
in the density of states and the absorption spectrum, it 
is clear that it plays a crucial role in the low-temperature 
dynamics of the excitons, because the dynamics is gov- 
erned by the local structure of wave functions and spec- 
tral distribution. 

As is observed in Fig. [5J above the states that extend 
over a few adjacent segments, states occur that are ex- 
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tended over many segments; they hardly carry any oscil- 
lator strength at all. In spite of this, these states play 
an important role in the problem of the temperature de- 
pendence of the exciton fluorescence decay time. Increas- 
ing the temperature, leads to their thermal population, 
which in turn leads to slowing down the fluorescence de- 
cay. 




mean mean 



FIG. 3: Distributions showing the statistics of the dimension- 
less oscillator strength per state, F v = Q^„=i fvn) 2 , collected 
in two narrow energy intervals of width 8E = 0.01 J. For the 
left panel this energy interval was chosen to be centered at the 
maximum of the spectrum of the oscillator strength per state 
(see Fig.^b)), while for the right panel this interval was cen- 
tered at the maximum of the absorption band. In each panel, 
solid (dashed) lines correspond to a = 0.1J (a — 0.5J). For 
each disorder strength, the distributions were collected using 
5000 random realizations of the disorder on linear chains of 
N = 500 sites. To stress the invariant nature of the distri- 
butions, the oscillator strengths have been rescaled by the 
average value F moall of F v in the interval and for the disorder 
strength under consideration. It is clearly seen that the width 
of each of the distributions is of the order of its mean. 

To end this section, we stress that the segment size 
iV* has the meaning of a typical number and in prac- 
tice undergoes large fluctuations^ First, the actual lo- 
calization length of exciton states, as may for instance 
be assessed from the participation ratio, is energy de- 
pendent; it becomes smaller towards lower energy and 
even over the narrow region of the absorption band, this 
difference may be sizablei^LS Second, even when we fo- 
cus on a narrow region within the absorption band, the 
stochastic nature of the disordered system gives rise to 
large fluctuations in the localization size of the exciton 
wave functions. For the optical dynamics, it is important 
that these fluctuations are also reflected in fluctuations of 
the oscillator strength of the exciton states. To illustrate 
this, we plot in Fig. |3Ia) the simulated distribution of 
the oscillator strength per state in an energy interval of 
width SE = 0.01J centered at the energy where the av- 



erage oscillator strength per state has its maximum (cf. 
Fig-IHb)). In Fig. E3[b), this is repeated with the nar- 
row interval centered at the energy where the absorption 
spectrum has its maximum (cf. Fig.^c)). The solid lines 
give the distributions for a = 0.1 J, while the dashed lines 
correspond to a — 0.5 J. The oscillator strength is given 
in units of its mean value F mean within the energy interval 
under consideration and at the disorder value considered 
[for a = 0.1J (0.5J), F mcan — 48 (13) at the absorption 
peak and F me an= 57 (15) at the peak of the oscillator 
strength per state]. We make the following observations: 
(i) The distributions have a width that is comparable to 
their mean, confirming the large fluctuations in segment 
size and extent and shape of the wave functions, (ii) 
In both narrow energy intervals, states are found with 
hardly any oscillator strength, which indicates that at 
both energies states occur with nodes. From the distri- 
butions it appears that such states with nodes are more 
abundant at the center of the absorption band (i.e., to- 
wards higher energy), which is consistent with the above 
explained local structure. We have confirmed this by 
doing a statistical analysis of the node structure of the 
wave function, using the value of | J2 n L P^n\ ( Pvn\\ (cf. the 
criterion for the node structure developed in Ref. I35T) . 
(iii) Interestingly, the distributions plotted on the scale 
of i^mean appear to be invariant under changing the disor- 
der strength. This is consistent with an invariance of the 
local structure of the exciton energies and wave functions 
near the lower band edge. 

III. INTRABAND RELAXATION MODEL 

To describe the dynamics of the exciton eigenstates \v) 
of the Hamiltonian with site disorder, we will account for 
the effect of spontaneous emission and the scattering of 
excitons on lattice vibations. We will describe this com- 
bined dynamics on the level of a Pauli master equation 
for the populations P v (t) of the exciton states: 

N 

Pu = -luPu + ^ W ^ P ^~ W ^ P -)^ ( 6 ) 

where the dot denotes the time derivative, 7„ = 
7o(X2^Li fun) 2 = JoFv is the spontaneous emission rate 
of the exciton state | v) (with F u the dimensionless oscil- 
lator strength), and W vji is the transition rate from the 
localized exciton state to the state \ v). This transition 
rate is crucial in the description of the redistribution of 
the exciton population over the localized exciton states. 
The model for W v>1 is based on certain assumptions about 
the coupling between the excitons and the lattice vibra- 
tions. In Refs. 13 and Il6l only the coupling between 
the excitons and the vibrations of the linear chain itself 
was taken into account. In reality, however, the exciton 
chains that we are interested in, such as linear aggregates, 
are not isolated, but are embedded in a host medium, so 
that the excitons are coupled to the host vibrations as 
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well. The density of states of the latter is large com- 
pared to that of the vibrations of the chain itself. For 
this reason, it is natural to assume that it is the coupling 
of the excitons to the host vibrations that determines the 
exciton intraband relaxation in the linear chain. 

In this paper, we adopt the glassy model for W vtl which 
we have introduced in our previous publication dealing 
with homogeneous molecular aggreates. 19 This model as- 
sumes a weak linear on-site coupling of the excitons to 
acoustic phonons of the host medium and ignores corre- 
lations in the displacements of the different sites on the 
chain. Within this model, the transition rates are given 
by (see for details Ref. Il9r) 

N 

W VI1 = Wo S{E U - Vln^ln 

n=l 

f n(E v -E^, E u >E, t , 

\ l + n(£! M -.E„), E v <E tl . 

Here, the constant Wq is a parameter that character- 
izes the overall strength of the phonon-assisted exciton 
scattering rates. Its microscopic definition involves the 
nearest-neighbor excitation transfer interaction J, the 
velocity of sound in the host medium, and the mass of 
the sites in the exciton chain^S In this paper, we will 
consider Wq as one composite parameter, which may 
be varied to account for different host and (or) exci- 
ton systems. The spectral factor S{E V — E^) describes 
that part of the \E V — E^-dependence of W Vfl which de- 
rives from the exciton-phonon coupling and the density 
of phonon states. The sum over sites in Eq. J7J) repre- 
sents the overlap integral of exciton probabilities for the 
states \n) and \v). Finally, n(fl) — [exp(f2/T) — is 
the mean occupation number of a phonon state with en- 
ergy O (the Boltzmann constant is set to unity). Due 
to the presence of the factors n(f2) and 1 + n(O), the 
transition rates meet the principle of detailed balance: 
W V n = W^v exp^-E^ - E v )/T\. Thus, in the absence of 
radiative decay (p/ v = 0), the eventual (equilibrium) ex- 
citon distribution is the Boltzmann distribution. 

Within the Debye model for the density of phonon 
states, the spectral factor is given byi^ S{E V — E^) = 
{\E V — E^\/J) 3 . It is worth noting, however, that the ap- 
plicability of this model to glassy host media is restricted 
to a very narrow frequency interval of the order of several 
wave numbers (see, for instance, Refs.l39land l40|) . There- 
fore, we rather consider a simple linear approximation for 
this factor, S(E L/ — E„) = \E V — E,,\l J, which is similar 
to the dependence used in Refs. [13 and yfl- This scaling 
properly accounts for a decrease of the exciton-phonon 
interaction in the long-wave acoustic limi'b^i^ and pre- 
vents the divergence of W vll at small values of \E V — E^ \ . 
We checked that the results reported in this paper are 
not essentially affected by assuming a higher power in 
the dependence of W vtl on the energy mismatch. 

In disordered systems, the overlap integral 
J2n=i 'filnffm appearing in the expression for W vfl 
plays a much stronger role in the optical dynamics 



than the details of the dependence of S on the energy 
mismatch. The fact that W vll is proportional to this 
overlap integral allows one to distinguish between two 
types of exciton transitions occurring in the vicinity 
of the bottom of the exciton band (the region of main 
interest at low temperatures), namely intrasegment 
and intersegment transitions. These two types of 
transitions involve, respectively, states localized on 
the same localization segment of the chain and states 
localized on different segments. As has been established 
in Ref. |3(| for both processes the overlap integrals scale 
inversely proportional to the typical segment size N*, 
while numerically the intrasegment overlap integral is 
approximately 50 times as large as the intersegment one 
(independent of the disorder strength a). Furthermore, 
the overlap integrals between the local states of a 
segment and one of the higher states that is extended 
over the same segment as well as the adjacent ones (such 
as the states 12, 13, 14, and 15 in Fig.|21), are of the same 
order as the intrasegment overlap integral. This implies 
a specific scenario for the exciton intraband relaxation 
at low temperatures. Let us assume that the exciton, 
initially created in the blue tail of the absorption band (a 
condition which is usually met in experiments), quickly 
relaxes to one of the exciton states of a local manifold. 
Denote the two local states involved by the quantum 
labels 1 (ground state) and 2 (excited state). Because 
of the difference in intra- and intersegment overlap 
integrals discussed above, the exciton first relaxes within 
the local manifold, provided that the intrasegment 
transition rate is larger than the radiative rates 7 M 
(/j = 1,2). Only after this first step of relaxation, the 
exciton can hop to the states of other (adjacent) local 
manifolds, again, provided the intersegment hopping 
rate is larger than 7^. Note that if 7 M is larger than any 
transition rate W vlll the exciton emits a photon before 
any relaxation step (intra- or intersegment) occurs. On 
this basis, one may distinguish between a fast and a 
slow limit of intrasegment relaxation. To this end, we 
compare the typical value of the intrasegment transition 
rate Wrz at zero temperature, with the radiative rate 
7i of the local ground state. If W\2 > 71, we are in the 
fast-relaxation limit. This inequality guarantees that 
at zero temperature, the exciton fluorescence decay is 
governed by the radiative process with a rate of the 
order of 71. By contrast, if W\% < 71, the intrasegment 
relaxation rate W12 dictates the exciton fluorescence 
decay. This is very similar to the distinction of fast and 
slow relaxation which we have previously made in the 
case of fully delocalized excitons in homogeneous molec- 
ular aggregates piS The reader may find an elaborate 
discussion on the principle difference between these two 
limits in Ref. 

For the remainder of this paper, it is useful to find 
the value of Wq that distinguishes between the limits 
of fast and slow intrasegment relaxation; we will denote 
this value by W(5 ntra . We first estimate W12 by replacing 
En^nfln ^ W, giving W 12 = W (E 2 - E 1 )/JN*. 
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Then we equate W12 to the superradiant decay rate 71, 
which is typically JqN*, and obtain 

= ■ (7) 

If we restrict for a moment to nearest-neighbor inter- 
actions, we may even go one step further by realizing 
that, as a consequence of wave vector quantization and 
level repulsion within localization segments, the typical 
energy difference E2 — E\ then scales like E2 — E\ w 
37r 2 J/jV* 2 t 17 i 35 i 43 which leads to 

W^ ntra = (Stt 2 )-^^* 4 . (8) 

From this expression, it is clear that W(5 ntra steeply de- 
creases with decreasing segment size, i.e., with increasing 
degree of disorder. This steep scaling is not essentially 
affected by taking into account long-range dipole-dipole 
interactions. Thus, if the degree of disorder increases 
and the excitons become more localized, there is a strong 
tendency of the system to move into the fast-relaxation 
limit, where the radiative decay governs the fluorescence 
kinetics. 

Following intraband relaxation, the next relaxation 
step involves transitions between states of different seg- 
ments (exciton migration). Similarly to the above, one 
may distinguish the limits of fast and slow intersegment 
relaxation, defined as Wy\ > 71 and Wy\ < 71, re- 
spectively. Here, the states I and 1' denote the ground 
states localized on adjacent segments and it is assumed 
that Ei > Ey . If the transfer process is slow, the ex- 
citon spontaneously decays before it makes a hop to an 
adjacent localization segment with Ey < E\. As a re- 
sult, the fluorescence spectrum is expected to coincide 
with the absorption band, because the latter is mainly 
determined by the total collection of local ground states 
(the segment states with appreciable oscillator strength) , 
i.e., in the absorption spectrum each local ground state 
contributes, independent of the value of its energy. How- 
ever, if the intersegment transfer process is fast, the ex- 
citon can make down-hill hops before the fluorescence 
is emitted. Therefore, lower-energy local ground states 
will give a larger contribution to the fluorescence spec- 
trum than higher-energy ones; this gives rise to a visi- 
ble red shift (Stokes shift) of the fluorescence spectrum 
with respect to the absorption band (see Sec. lIVB|) . Us- 
ing the above noted factor of 50 difference between the 
intra- and intersegment overlap integrals^ the value of 
Wq that distinguishes between the limits of fast and slow 
intersegment relaxation, denoted by Wg nter , is given by 
W inteI = 5(W,j ntra . With reference to Eq. ©, we note 
that, keeping Wo constant, a Stokes shift is expected to 
become more noticeable for smaller segments, i.e., for 
growing disorder strength. 



IV. STEADY-STATE FLUORESCENCE 
SPECTRUM 

In this section, we deal with the exciton dynam- 
ics and the corresponding fluorescence spectrum under 
steady-state conditions, which are maintained by opti- 
cally pumping the system. The steady-state fluorescence 
spectrum is defined as 




where Pjf is the solution of the steady-state master equa- 
tion 

N 

= R V - lv Pf + ^(W vtt If - W^Rf) ■ (10) 

Here, R v denotes the constant rate of optically creating 
population in the vtYi exciton state by a c.w. pump pulse. 

We performed numerical simulations of exciton fluo- 
rescence spectra, by solving the master equation Eq. I|I(J|) 
for randomly generated realizations of the disorder. For 
each realization, we diagonalize the Hamiltonian H nm 
in order to calculate the radiative constants 7^ and the 
transition rates W v)1 entering the master equation. In all 
simulations reported in the remainder of this paper, we 
set J = 600 cm" 1 (1.8 x I0 13 s" 1 ) and 7o = 2 x I0~ 5 J 
(3.6 x I0 8 s _1 ). These parameters are quite typical for 
molecular aggregates of polymethine dyes, such as PIC. 

A. Steady-state fluorescence spectra at zero 
temperature 

In this subsection, we concentrate on the zero- 
temperature fluorescence spectra, which may be nicely 
used to illustrate the different regimes of exciton relax- 
ation. In Fig. 0] we present by dashed and solid lines 
the zero-temperature steady-state fluorescence spectra 
for various combinations of values of the disorder strength 
a and the phonon-assisted exciton scattering strength 
Wo- All spectra were calculated under the condition of 
off-resonance optical pumping in a narrow energy window 
of width 0.05 J in the blue tail of the absorption band. 
The exact position of the pump window was chosen to 
be blue-shifted relative to the maximum of the absorp- 
tion band (simulated at the same value of the disorder 
strength) by three times the full width at half maximum 
(FWHM) of this band. The pump rate of each exciton 
state v inside the pump window was taken proportional 
to its oscillator strength: R v = F v . As is clear from 
Fig. 0] for each value of a, the fluorescence spectrum 
nearly follows the absorption band (shown by the dotted 
line) as long as Wq is small, while this spectrum expe- 
riences a visible Stokes shift if Wq is increased. This 
agrees with our expectations formulated in Sec. lIIII Also 
in agreement with our arguments made at the very end 
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of Sec. IIIII we observe that (for constant Wo) the mag- 
nitude of the Stokes shift is smaller for smaller disorder 
strength. 
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FIG. 4: Zero-temperature steady-state exciton fluorescence 
spectra calculated for various disorder strengths a and exciton 
scattering rates Wq. Spectra were obtained by numerical so- 
lution of the steady-state master equation Eq. 1101 under the 
condition of off-resonance blue-tail optical pumping (see text 
for details), (a) a = 0.1J, with W = J (dashed), Wo = 100J 
(dash-dotted), and Wo = 10 5 J (solid); (b) a = 0.3 J, with 
Wo = 0.1J (dashed) and W = 100 J (solid); (c) a = 0.5 J, 
with Wo = 0.01J (dashed) and Wo = 100 J (solid). The dot- 
ted line in each panel represents the absorption band, while 
the solid vertical line shows the location of the maximum of 
the spectrum for the oscillator strength per state. The other 
parameters used in the simulations were N = 500, J = 600 
cm -1 , and 70 = 2 x 10 -5 J. The average was performed over 
5000 realizations of the disorder, using energy bins of 0.005J 
to collect the fluorescence spectrum. 

To gain more quantitative insight into the above be- 
havior of the fluorescence spectra, let us estimate the 
values of the parameters Wo ntra and Wq 11 * 61 that distin- 
guish the limits of fast and slow relaxation for intra- and 
intersegment transitions, respectively, as introduced in 
Sec. lIIII First, we do this for the largest disorder strength 
considered in the simulations, a = 0.5 J (Fig.QJc)). From 
the maximum of the spectrum for the oscillator strength 
per state plotted in Fig. IHb), we find as typical seg- 
ment size N* = 15 for this disorder strength. Similarly, 
the typical separation between the two bottom states 
of a localization segment may be estimated from the 
FWHM of the absorption band, giving E 2 — Ei = 0.4J 
(Fig. nfc)). Substituting these data into Eq. Q, we 
obtain W<j ntra « 0.01J and W inter = 50Wf] ntra » 0.5J. 
Thus, the smaller value of Wo considered in the simula- 
tions (0.01 J) is equal to W lIrtra , while it is much smaller 
than Wo nter , i.e., no intersegment hops will occur prior to 
fluorescent emission. The latter explains why there is no 
Stokes shift of the fluorescence spectrum. On the other 
hand, as the system is in the intermediate regime with re- 
gards to the intrasegment hopping (Wo = Wo" tra ), the ra- 
diative channel can compete with the intrasegment relax- 
ation. This explains the presence of the small and narrow 
fluorescence peak coinciding with the pumping interval 
in the blue tail of the absorption spectrum. By contrast, 
the higher value of W = 100J exceeds W intcl = 0.5J 
by more than two orders of magnitude. This results in a 
visible Stokes shift of the fluorescence spectrum as well 
as a strong reduction of the fluorescence in the excitation 
window. 



Analogous estimates performed for a = 0.1J 
(Fig. Ufa)) bring us to the following results: WQ lltla ss J 
and W,j nter = 50W,j ntra « 50J. Here we used N* = 57 
and E2—E1 = 0.04 J, taken as previously from Figs.^b) 
and (c), respectively. For the smaller value of Wo = J, 
the excitons are again within the intermediate regime 
with regards to the intrasegment relaxation. As a con- 
sequence, the fluorescence spectrum shows the same pe- 
culiarities as in the case of a — 0.5J at Wq = 0.01J. 
However, as the higher value of Wo = 100 J is only twice 
as large as compared to W lnter ~ 50J, the Stokes shift 
here is smaller than in in the case of the higher disorder 
magnitude. A large Stokes shift may be forced by taking 
W = 10 5 J, as is also illustrated in Fig. Ufa). It should 
be noted, however, that this large value for Wo lies out- 
side the range of validity of our theory. The reason is 
that such a large scattering rate leads to W12 3> E% — E\ 
(W12 ~ 0.8J for the current example), which implies 
that the second-order treatment of the exciton-phonon 
interaction is a poor approximation. More importantly, 
under such conditions the segment picture breaks down, 
because the exciton coherence size is dominated by the 
scattering on phonons instead of static disorder. A 
proper description then requires using a density matrix 
approach 44 ! 45 



B. Temperature dependence of the Stokes shift 

We now turn to the temperature dependence of the 
steady-state fluorescence spectrum. In particular, we are 
interested in the temperature dependence of its Stokes 
shift with respect to the absorption band. We have calcu- 
lated this shift as the difference in peak positions between 
the absorption and fluorescence bands. As the simu- 
lated fluorescence spectrum contains appreciable stochas- 
tic noise (it is not possible to apply the same smoothen- 
ing as may be used when simulating the absorption 
spectrums») , its peak position was determined by fitting 
the upper half of the peak to a Gaussian lineshape. Fig- 
ure shows the thus obtained results, for three differ- 
ent values of Wo at a fixed disorder strength, a = 0.3 J. 
The characteristic peculiarity of all curves is that they 
are non-monotonic: the Stokes shift first increases upon 
heating and then goes down again. The extent of the 
temperature interval over which the Stokes shift increases 
is small compared to the absorption band width, which 
is 170 K for the current degree of disorder. A Stokes 
shift that increases with temperature is counter-intuitive, 
because, at first glance, it seems that the temperature 
should force the excitons to go up in energy, giving rise to 
a monotonic decrease of the Stokes shift. This expected 
behavior is indeed observed in inhomogeneously broad- 
ened systems doped with point centers like, for instance, 
glasses doped with rare-earth ions4& The explanation for 
the peculiar behavior of the Stokes shift in the Frenkel 
exciton chain is similar to that for the nonmonotonic be- 
havior found for disordered quantum wells ^ and resides 
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in a thermally activated escape from local potential min- 
ima. However, for the case of the disordered linear chain, 
the detailed knowledge of the low-energy spectral struc- 
ture provides additional means to unravel the character- 
istics of this behavior. 
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FIG. 5: Semi-log plots of the temperature dependence of 
the Stokes shift of the fluorescence spectrum at the disor- 
der strength a = 0.3J, for different exciton scattering rates: 
Wo = J (solid), Wo = 10 J (dashed), and W = 100 J (dot- 
ted). The data were obtained by numerical solution of the 
master equation Eq. HUH under the condition of off-resonance 
blue-tail excitation (see text for details). The dots mark the 
numerical data, while the lines connecting the dots are guides 
to the eye. The other parameters used in the simulations were 
N = 500, J = 600 cm" 1 , and 70 = 2 x 10~ 5 J. The average 
was performed over 5000 realizations of the disorder. 



Let us first recall the zero-temperature scenario of the 
exciton relaxation. Excitons created initially at the blue 
tail, rapidly relax to the local states of the DOS tail, 
which are visible in fluorescence. After that, the exci- 
tons may relax further, whithin the manifold of the local 
ground states. At zero temperature, however, this pos- 
sibility is very restricted. The reason is that an exciton 
that relaxed into one of the local ground states may move 
to another similar state of an adjacent localization seg- 
ment only when the latter has an energy lower than the 
former and W exceeds W™ tcr (Sec. llllf) . The typical en- 
ergy difference between the local ground states is of the 
order of the absorption band width. Therefore, after one 
jump the exciton typically resides in the red tail of this 
band. The number of states with still lower energy then 
strongly reduces, giving rise to an increased expectation 
value for the distance to such lower energy states. In 
fact, already after one jump the exciton has a strongly 
suppressed chance to jump further during its lifetime; it 
will generally emit a photon without further jumps (mi- 
gration). Thus, the states deep in the tail of the DOS can 
typically not be reached by the excitons, simply because 
they occur at a low density. This qualitatively explains 
the fact why the Stokes shift of the fluorescence spec- 
trum does not exceed the absorption band width, even 



for Wo large compared to WQ ltcr (see Fig. |2J. Upon a 
small increase of the temperature from zero, however, it 
becomes easier to reach those lower-lying states, because 
the spatial migration to other segments may take place by 
thermally-activated transitions involving exciton states 
that are extended over several localization segments as 
intermediate states4£ It is this indirect hopping that is 
responsible for the increase of the Stokes shift at tem- 
peratures small compared to the absorption band width. 
Further heating will thermalizc the excitons and lead to 
real populations of higher-energy states; the Stokes shift 
will then decrease again. 

To be more specific, we present estimates. We first 
note that the overlap integral of the squared wave func- 
tions in Eq. J7J) for a local ground state 1 and a higher 
state 3 that extends over more than one segment, but 
still overlaps with the ground state 1, has the same order 
of magnitude as for the states f and 2 of a local manifold, 
i.e., 1/N* . Then, the transition rate up from the ground 
state to the more extended state, W31, can be estimated 
as W31 « [W (E 3 - S 1 )/J7V*]exp[-(S 3 - E{)/T]. In 
order for the exciton migration via the higher state to 
be activated, W31 should be larger than 71 = 70-/V* , the 
spontaneous emission rate of state 1. Equating these two 
rates gives us a temperature To at which the Stokes shift 
is increased over its zero-temperature value: 

. loN* 2 J 

Taking as an estimate for E 3 — E\ the FWHM of the 
absorption band, 0.2 J at a = 0.3 J, we obtain T = 57 K, 
32 K, and 22 K for W = J, 10J, and 100J, respectively. 
These numbers are in good agreement with the positions 
of the maxima of the curves in Fig. [S] 

To the best of our knowledge, the only one-dimensional 
exciton system for which the temperature dependence of 
the Stokes shift has been measured, is the molecular ag- 
gregate that is formed by the cyanine dye THIATS. 15 
For this aggregate, indeed a non-monotonic temperature 
dependence of the Stokes shift was reported at low tem- 
peratures, very similar to our numerical results. Thus, 
the model we are dealing with provides an explanation 
of the behavior reported in Ref. ^[ A detailed fit to 
these experimental data, including also the absorption 
spectrum and the fluorescence lifetime of this aggregate, 
will be presented in a separate paper4& 



V. FLUORESCENCE DECAY TIME 

We proceed to study the decay time of the total time- 
dependent fluorescence intensity following pulsed excita- 
tion at t — 0. This intensity reads 



T = (E 3 -E 1 )/]n 



(11) 



I(t) 



(12) 



where the P„ (t) are obtained from the Pauli master equa- 
tion Eq. © with the appropriate initial conditions. It 
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should be realized that the proper definition of the de- 
cay time requires attention, as the time dependence of 
I(t) is not mono-exponential. We already encountered 
this problem for homogeneous aggregates^ but for the 
disordered systems under consideration, the problem is 
even more obvious. The multi-exponential behavior is a 
consequence of the large fluctuations in the spontaneous 
decay rates of different exciton states, as is clearly demon- 
strated by the distribution of exciton oscillator strengths 
plotted in Fig. [3] Similarly, large fluctuations occur in the 
transition rates W„ M . The simplest solution is to define 
the decay time, r e , as the time it takes the fluorescence 
intensity I(t) to decay to 1/e of its peak value I(t pea k)'- 

I(tpeak + T e ) = -I(t pea k) ■ (13) 

Throughout this paper, we will use this definition of 
the decay time. It should be noted that in the case 
of initial excitation in the blue tail of the absorption 
band, t pea k 7^ 0, due to that fact that a finite time 
elapses before the exciton population reaches the lower- 
lying emitting states. An alternative choice for the de- 
cay time is the expectation value of the photon emission 
time, t = Jq 00 Pv(t)) dt. For mono-exponential de- 
cay, both definitions give the same result, but in general 
this does not hold. In particular, for non-exponential 
fluorescence kinetics, the latter definition, t, only gives a 
meaningful measure of the fluorescence time scale in the 
limit of the fast intrasegmcnt relaxation (sec for more 
details the discussion presented in Ref. ITfllh 



A. Broadband resonance excitation 

We first consider the case of broadband resonance ex- 
citation, which is similar to what takes place in echo ex- 
periments^ Under this condition, all states are excited 
with a probability that is proportional to their respec- 
tive oscillator strengths, P v [t = 0) = F V: meaning that 
the spectral profile of the initially excited states coin- 
cides with the absorption band. Thus, the initial exciton 
population mostly resides in the superradiant states. 

In Fig. EJa)-(c), we depicted the temperature depen- 
dence of the fluorescence decay time r e for three different 
disorder strengths, a — 0.1J, 0.3J, and 0.5J, respectively; 
for each a value two different strengths of the exciton 
scattering strength Wo were considered. The solid line 
in each panel presents results for the intermediate regime 
with regards to the intrasegment relaxation, i.e., when 
W21 ~ 71, while the dashed line shows results in the 
limit of fast relaxation, W21 3> 71 . For all three disorder 
strengths, the higher value of Wq considered, is below the 
threshold W~Q nteI for fast intersegment relaxation, in other 
words, no visible Stokes shift occurs in the fluorescence 
spectra for any of the chosen parameters. The insert in 
Fig. Eta) shows on a semi-log scale the time-dependence 
of the fluorescent traces underlying the reported decay 
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FIG. 6: Temperature dependence of the fluorescence decay 
time r e calculated for various disorder strengths a and exciton 
scattering rates Wo- The data were obtained by numerical 
solution of the master equation Eq. @ under the condition of 
broadband resonant excitation, setting P v {t = 0) = F v . The 
dots mark the numerical data, while the lines connecting the 
dots are guides to the eye. (a) a — 0.1J, with Wo = J (solid) 
and Wo = 10 J (dashed); (b) a = 0.3 J, with W = 0.1J 
(solid) and W = J (dashed); (c) a = 0.5 J, with Wo = 0.01 J 
(solid) and Wo = J (dashed). The other parameters used in 
the simulations were N = 500, J = 600 cm" 1 , and 70 = 
2 x 10~ 5 J. The average was performed over 50 realizations 
of the disorder. The insert in (a) shows the time-dependence 
of the fluorescence intensity for Wo = 10J at four different 
temperatures: from top to bottom, the curves correspond to 
T = K, 17 K, 34 K, and 84 K, respectively. 



times for Wq = 10 J. These traces clearly show that in 
general the intensity decay is non-exponential. 

It is worthwhile to estimate the zero-temperature val- 
ues of the fluorescence decay time using the relationship 
r e = 1/(70 JV*). As previously, we take for N* the maxi- 
mum value of the spectrum for the oscillator strength per 
state (Fig. Eft))), i.e., N* = 57, 23, and 15 for a = 0.1 J, 
0.3J, and 0.5J, respectively. Then, for the parameters 
used in our simulations, 70 = 2 x 10~ 5 J and J = 600 
cm" 1 , the corresponding values of r e are 49, 121, and 
185 ps, respectively. As is seen, these estimated times 
are larger than the calculated ones in Fig. El The rea- 
son for this deviation is the resonance excitation condi- 
tion, combined with the large fluctuations in the oscilla- 
tor strength per state (Fig. EJ). Indeed, the states with 
a higher than typical oscillator strength are excited to a 
larger extent than those with typical (and smaller) os- 
cillator strengths. Obtaining a relatively large part of 
the initial population, they mainly determine the initial 
stage of the fluorescence kinetics. This gives rise to a 
faster decay rate than the typical one. 

Apart from some low-temperature peculiarities of T e 
for higher degree of disorder, all curves in Fig. El tend 
to go down upon increasing the temperature. This 
has the following explanation. The up-hill transition 
processes, which are characterized by the rate W21 oc 
exp[— (E2— Ei)/T], come into play when the temperature 
is increased. At some disorder dependent temperature, 
W21 becomes larger than 71, and the exciton popula- 
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tion from the initially populated superradiant states is 
transferred to higher (dark, initially not excited) states. 
This nonradiative loss of population from the superradi- 
ant states gives rise to a drop in the fluorescence intensity, 
which contributes to the observed fluorescence decay. We 
stress that it is the rate W21 that determines the time of 
this uphill process (cf. Ref. Il9l) . Thus, with increasing 
temperature, the drop in the fluorescence decay time T e 
reflects in fact the shortening of the intrasegment relax- 
ation time and not the exciton radiative lifetime. 

To conclude this subsection, we note that the range of 
variation of the fluorescence decay time with temperature 
differs dramatically for different disorder strengths. In 
particular, for a = 0.1J, r e decreases from its maximal 
value of 36 ps (at T = 0) to about 1 ps at 150 K, while for 
a = 0.5J this drop consists of only 15-20% of the zero- 
temperature value of T e . This simply results from the 
fact that the absorption band widths for these two values 
of the disorder strength are (in temperature units) 40 K 
and 300 K, respectively. We recall that the rate of up-hill 
transfer W21 oc exp[-(£ , 2 — E\)/T], E 2 — E\ being of the 
order of the absorption bandwidth. Thus, for a = 0.1 J, 
a temperature of 40 K is already sufficient to activate 
the up-hill transfer of population and to noticeably drop 
the fluorescence intensity. By contrast, even the highest 
temperature considered in the simulations, T — 160 K, 
is not enough to start the up-hill process for a = 0.5 J. 



B. Off-resonance blue-tail excitation 

We now turn to the case of off-resonance excitation 
in the blue tail of the absorption band. This is the 
usual situation in fluorescence experimentsi£*£*£*i2*ii*i2*i^ 
We recall that in this case between the absorption and 
emission events an additional step exists: the vibration- 
assisted relaxation from the initially excited states to the 
radiating ones. This results in different scenarios for the 
exciton fluorescence kinetics, dependent on the relation- 
ship between the intraband relaxation rate and the rate 
of exciton spontaneous emission. 19 

In Fig. we depict the temperature dependence of the 
fluorescence decay time r e obtained from numerical sim- 
ulations for various strengths of the disorder a and the 
vibration- assisted exciton scattering rate Wq. The initial 
condition for solving Eq. © was taken P v (t = 0) = F v 
in a narrow window in the blue tail of the absorption 
band, defined in the same way as in Sec. If V Al The scat- 
tering rates Wq were chosen to realize different limits of 
the intraband relaxation. In particular, in the case of the 
smallest disorder strength, a = 0.1J, the scattering rates 
Wq = 0.1 J and Wo — J describe the limits of interme- 
diate and fast intrasegment relaxation, respectively (see 
the discussion presented in Sec. IIV Afl . However, with 
respect to the intersegment hopping, these values both 
correspond to the slow limit. Finally, the highest value 
of Wq — 100J describes the limit of fast intersegment 
relaxation. Similar relationships exist between the Wq 
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FIG. 7: Temperature dependence of the fluorescence decay 
time Te calculated for various disorder strengths a and exciton 
scattering rates Wo. The data were obtained by numerical so- 
lution of Eq. (|SJ under the condition of off-resonance blue-tail 
excitation (see text for details). The dots mark the numerical 
data, while the lines connecting the dots are guides to the eye. 
(a) a = 0.1J, with Wo = J (dotted), Wo = 10J (dashed), 
and Wo = 100 J (solid), (b) a = 0.3 J, with Wo = 0.1 J (dot- 
ted), Wo = J (dashed), and Wo = 10J (solid), (c) a = 0.5J, 
with Wo = 0.01J (dotted), W = 0.1J (dashed), and Wo = J 
(solid). The other parameters used in the simulations were 
N = 500, J = 600 cm" 1 , and 70 = 2 x 10 _5 J. The average 
was performed over 50 realizations of the disorder. 



values for the other degrees of disorder. 

From Fig. we observe that, in contrast to the case of 
the broadband resonance excitation, all r e curves go up 
almost linearly with temperature, independent of the val- 
ues for a and Wq ■ Being well separated at zero tempera- 
ture, they tend to approach each other at higher temper- 
atures. The latter effect is more pronounced for smaller 
degrees of disorder. Some of the curves show a low- 
temperature plateau, whose extent is smaller than the 
absorption bandwidth. On the other hand, there also is 
a common feature between the resonant and off-resonant 
type of excitation: the range of variation of r e with tem- 
perature is smaller for more disordered systems. This 
effect has the same explanation as discussed in Sec. IV Al 

Two processes are responsible for the observed increase 
of r e with temperature: intraband down-hill relaxation 
after the excitation event and thermalization of the exci- 
tons over the band. Let us consider the first step of the 
population transfer from the initially excited states to 
lower states, both dark and superradiant. An important 
feature of this process is that it is almost nonselective 
due to the linear dependence of the transition rates W Vfi 
on the energy mismatch. As a result, the lower states 
are populated almost equally, whether superradiant or 
dark. This is in sharp contrast to the case of resonance 
excitation where, in fact, only the superradiant states are 
initially excited. The further scenario can be understood 
by considering a simple two-level model. Let level 1 de- 
note the lowest, superradiant, state of the local manifold, 
having an emission rate 71, while level 2 is the higher- 
energy, dark, local state. We assume that following the 
fast initial relaxation described above, both levels are 
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excited equally. The Pauli master equation Eq. JH| now 
reduces to 

A = -(71 + W 2 i)Pi + W 12 P 2 ., (14a) 



P 2 = -W 12 P 2 + W21P1 , (14b) 

with initial conditions -Pi(O) = P 2 (0) — 1/2. We seek the 
solution of Eqs. HU) in the limits Wi 2 = W21 = W (T > 
E 2 — Ei) and W 71 (fast intrasegment relaxation). It 
is easy to find that the intensity / = —Pi — P 2 = 71 Pi is 
given by 

The second term in Eq. l|15f) can be neglected. From the 
first one it follows that the fluorescence decay rate is given 
by 7i/2, which directly reflects the exciton's radiative de- 
cay. It is only half the superradiant rate 71 due to the fast 
exchange of population between the superradiant level 1 
and the dark level 2. If the temperature is increased, 
one should generalize this discussion by considering the 
situation where I non-radiating levels, equally populated 
initially, are rapidly exchanging population with the su- 
perradiant level. The number I increases with tempera- 
ture. The result is straightforward: one should replace 
the rate 71/2 by 7i/(Z + 1). This qualitatively explains 
the temperature behavior of the exciton fluorescence de- 
cay time found in the simulations. 

To conclude the discussion of the numerical results pre- 
sented in Fig. we comment on the zero-temperature 
value of T e . For each a value, this noticeably depends 
on Wq, decreasing as Wo goes up. Furthermore, the cal- 
culated zero-temperature values for r e at the smallest 
magnitudes of Wo considered in the simulations (Wo = 
J, 0.1J, and 0.01J for a — 0.1 J, 0.3J, and 0.5J, respec- 
tively), are larger than those estimated from the maxi- 
mum of the spectrum for the oscillator strength per state 
(Fig. IHb)). Recall that these estimates are 49 ps for 
a = 0.1 J, 121 ps for a = 0.3J, and 185 ps for a = 0.5J 
(see Sec. EH - 

The observed decrease of the zero-temperature value 
of r e with Wo, may be understood from the fact that 
the emitting exciton sees a distribution of the oscilla- 
tor strength that differs for different Wo values. Below, 
we provide a qualitative picture of this. The exciton is 
initially excited at the blue tail of the absorption band, 
where the oscillator strengths are small, so that down-hill 
relaxation dominates over the emission. This allows the 
exciton to go down in energy until the intraband relax- 
ation rate becomes comparable to or smaller than the ra- 
diative rate. Once this has happened the exciton emits a 
photon. Let us analyze first what happens at the smallest 
magnitudes of the exciton scattering rate Wq considered 
in the simulations (see above). These values correspond 
to the intermediate case with regards to intrasegment re- 
laxation, while with respect to the intersegment hopping, 



they fall in the slow limit. This means that the exciton, 
after it has relaxed from the blue-tail states to the su- 
perradiant states, does not move any more. It can then 
only emit a photon. The zero-temperature steady-state 
spectra presented in Fig.0]providc information about the 
spectral location of the exciton emission. For the values 
of Wo we are discussing, the emission spectra coincide 
in general with the absorption spectra. At the same 
time, the maximum of the oscillator strength distribu- 
tion is shifted to the red from the absorption maximum 
(cf. Fig.^). Therefore, the fluorescence decay time is ex- 
pected to be larger than that estimated via the maximum 
of the oscillator strength distribution. This explains the 
results of the simulations. 

For the largest value of Wo considered in the simula- 
tions (W = 100J, 10J, and J for cr = 0.1J, 0.3J, and 
0.5 J, respectively), the limit of fast intersegment relax- 
ation applies. This means that after the fast intraseg- 
ment relaxation to the superradiant states, the exciton 
still has a chance to relax further due to intersegment 
hops (migration). As a result, the emission spectra are 
shifted towards the maxima of the spectra of the oscil- 
lator strength per state (see Fig. @}, which explains the 
shortening of the decay time observed in the numerical 
simulations with increasing value of Wo. 

C. Dependence on the detection energy 

To end our analysis of the fluorescence decay time, we 
address its dependence on the detection energy. This has 
attracted considerable attention in the literature on ag- 
gregates and polymers (see for instance Refs. I15TI23T) . To 
study this dependence, we have simulated the detection 
dependent fluorescence intensity, defined through: 

I(E d ;t) = (^2j u P u (t)A(E d - E u )j , (16) 

where Ed denotes the central detection energy and 
A(Ed — EJ) is the detection window, which is unity for 
\Ed — E v \ < 0.0025 J and zero otherwise. Simulations 
were carried out for aggregates of N = 250 molecules 
with cr = 0.3J and Wo = 100J, other parameters as 
usual. We have considered blue-tail short-pulse excita- 
tion conditions, as was done in Sec. IV Bl and three de- 
tection energies: the peak of the steady-state fluorescence 
spectrum (cf. Fig. Hfb)), and the positions of the blue 
and red half maximums of this spectrum. From the in- 
tensity traces, we have extracted the energy dependent 
1/e decay times T e (Ed). The results as a function of tem- 
perature are shown in Fig. [HJ with the upper, middle, and 
lower curve corresponding to red, peak, and blue detec- 
tion energy, respectively. 

We see that at zero temperature the decay times clearly 
differ for the three detection energies. For the case of red 
detection, we observe a decay time of 152 ps. Based 
on the average oscillator strength per state at this red- 
wing energy, we arrive at a purely radiative decay time 
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FIG. 8: Temperature dependence of the fluorescence decay 
time T e {Ed) calculated for aggregates of N = 250 molecules, 
with J = 600 cm" 1 , 70 = 2 x 10" 5 J, a = 0.3J, and W = 
100 J. The data were obtained by numerical solution of Eq. © 
under the condition of off-resonance blue-tail excitation (see 
text for details). The dots mark the numerical data, while 
the lines connecting the dots are guides to the eye. The three 
curves correspond to different detection energies Ed related to 
the steady-state fluorescence spectrum in Fig.0Jb): detection 
at the position of the red half maximum of this spectum (solid 
line), at the peak of this spectrum (dashed line), and at the 
blue half maximum (dotted line). The average was performed 
over 10000 realizations of the disorder. 



(l/ r yoF(E)) of 145 ps. The agreement between these 
numbers indicates that at the red side the decay time at 
zero temperature is determined completely by radiative 
decay; intraband relaxation has no effect at this energy. 
The reason is that this detection energy lies very deep 
in the tail of the DOS, where the occurrence of neigh- 
boring segments with lower energy is negligible. For the 
detection at the peak position, we find a zero temper- 
ature decay time of 100 ps, which is about 16% faster 
than the purely radiative time scale of 120 ps at this en- 
ergy. The difference is due to relaxation to lower lying 
exciton states in neighboring segments. Naturally, this 
effect is even stronger at the blue position, where the 
observed decay time of 44 ps is considerably faster than 
the purely radiative decay time of 168 ps. At this blue 
position, also an increased influence from intrasegment 
relaxation exists, as at these higher energies a fraction 
of the states already represent an excited segment state, 
with one node (cf. Fig.^J). 

We also observe from Fig. |H1 that at high temperature 
the decay times for the different detection energies ap- 
proach each other and in fact they then all tend to the de- 
cay time of the total fluorescence intensity (cf. Fig.[J£b)). 
This is a consequence of the fact that the scattering rates 
are then large enough for the exciton populations to be- 
come equilibrated on the time scale of emission. A com- 
parable observation was made in Ref. for the decay 
of the populations of the various exciton states for ho- 
mogeneous aggregates. We finally observe that the tem- 



perature dependence is non-monotonic in the same tem- 
perature range for which the Stokes shift behaves non- 
monotonically (cf. Fig. 0. Indeed, we attribute this be- 
havior to the same temperature activated intersegment 
relaxation via higher lying exciton states. At the blue 
detection side (which still lies in the red tail of the DOS), 
this effect leads to a decrease of the lifetime, as it opens 
extra decay channels. On the red side, the situation is 
more subtle. This energy is so deep in the red tail of the 
DOS, that even the activated migration hardly opens new 
channels for decay. Instead, the activated relaxation oc- 
curring at the blue side towards lower energies will be 
lead to extra contributions in the fluorescence intensity 
at the red side and, thus, to a growth of the decay time 
at this energy. At the peak position, we deal with the 
intermediate situation and we see a very small net effect. 

We notice that the general characteristics displayed in 
Fig. [S] very well cover the experimental results reported 
by Scheblykin et alJ£ 



VI. SUMMARY AND CONCLUDING 
REMARKS 

We have performed a numerical study of the tempera- 
ture dependence of the exciton dynamics in linear Frenkel 
exciton systems with uncorrelated diagonal disorder. In 
particular, we have focused on the resulting temperature 
dependent steady-state fluorescence spectrum, its Stokes 
shift relative to the absorption spectrum, and the decay 
time of the total fluorescence intensity following pulsed 
excitation. The complicated exciton dynamics reflected 
in these observables is governed by the interplay between 
thermal redistribution of the excitons over a set of eigen- 
states, which are localized by the disorder, and their ra- 
diative emission. The redistribution of exciton popula- 
tion within the manifolds of localized exciton states was 
described by a Pauli master equation. The transition 
between two localized states was assumed to originate 
from the coupling of the excitons to acoustic phonons of 
the host medium, the transition rates being proportional 
to the overlap integral of the corresponding wave func- 
tions squared. The model is characterized by two free 
parameters, er, which denotes the degree of disorder, and 
Wo, the phonon-assisted exciton scattering rate, which 
sets the overall scale for transition rates between exciton 
states. 

The fact that our model only accounts for scattering on 
acoustic phonons, in principle limits us to temperatures 
of the order of 100 K and less. This covers the tem- 
perature range of many experiments performed on linear 
dye aggregates. Besides, as is clear from our simulations, 
if the exciton scattering rate is large enough, equilibra- 
tion within the exciton space (on the time scale for ra- 
diative dynamics) already occurs at temperatures below 
100 K. Above the equilibration temperature, the precise 
nature of the scattering mechanism becomes unimpor- 
tant. In practice, the fact that in our Pauli master equa- 
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tion for the exciton populations, homogeneous broaden- 
ing (dephasing) can not be considered, probably yields 
a stronger limitation to the accuracy at elevated tem- 
peratures than the restriction to acoustic phonons. In- 
corporating dephasing, i.e., accounting for the possible 
breakdown of coherence within derealization segments, 
requires considering the exciton density operator (also 
see end of Sec. II V A)) . 

^,From our simulations we found that the Stokes shift 
of the fluorescence spectrum shows an anomalous (non- 
monotonic) temperature dependence: it first increases 
upon increasing the temperature from zero, before, at a 
certain temperature, it starts to show the usual mono- 
tonic decrease. This behavior was found previously for 
disordered quantum wells and physically derives from 
thermal escape from local potential minimal* We have 
shown that for disordered chains, the details of this 
behavior and the temperature range over which the 
anomaly takes place can be understood from the specific 
features of the exciton energy spectrum in the vicinity of 
the lower band edge: it is formed of manifolds of states 
localized on well separated segments of the chain and 
higher states that are extended over several segments. 
The migration of excitons between different segments 
augmented via intermediate jumps to higher states, is 
responsible for the non-monotonic behavior found in our 
simulations. Interestingly, such a non-monotonic behav- 
ior of the Stokes shift has recently been observed for the 
linear aggregates of the cyanine dye THIATS in a glassy 
host>i£ 

We also found that the temperature dependence of 
the decay time of the total fluorescence intensity is very 
sensitive to the initial excitation conditions. For broad- 
band resonant excitation, the fluorescence decay time de- 
creases upon increasing the temperature. The reason 
is that the initially created population of superradiant 
states is transferred to higher (dark) states. It is the 
time of this transfer that determines the fluorescence de- 
cay time. As this transfer time decreases with growing 



temperature, a decreasing fluorescence lifetime is found. 
Because fluorescence experiments are hard to perform 
under resonant excitation, it will be difficult to observe 
this effect of intraband redistribution in fluorescence. It 
would be of interest, however, to study its effect on res- 
onantly excited photon echo experiments. 49 

In the case of off-resonance blue-tail excitation (the 
condition that is usually met in fluorescence decay exper- 
iments), the fluorescence decay time goes up with grow- 
ing temperature, showing a nearly linear growth after a 
low-temperature plateau. The extent of the plateau de- 
pends on both the absorption bandwidth and on the ratio 
of the rates for exciton hopping and radiative emission. 
This behavior, a decay time that grows with tempera- 
ture, with a possible plateau at low temperatures, agrees 
with fluorescence experiments performed on the J-bands 
of linear molecular aggregates. It is of particular interest 
to note that a nearly linear dependence has been ob- 
served for aggregates of BIG— and THIATS 1 ^. Based on 
the density of states of homogeneous exciton systems, it 
has been suggested that such a linear dependence could 
only occur for two-dimensional systems^ It follows from 
our results that in the presence of disorder, a ubiquitous 
ingredient for molecular aggregates, one-dimensional ex- 
citon systems may exhibit such a linear temperature de- 
pendence as well. In a separate paper^ we will show 
that using the model analyzed in the present paper, it is 
possible to obtain good quantitative fits to the absorp- 
tion spectrum, the temperature dependent Stokes shift, 
as well as the temperature dependent fluorescence life- 
time measured for aggregates of the dye THIATS 
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